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TREE BASED FUNCTIONAL EXPANSIONS FOR FEYNMAN KAC 

PARTICLE MODELS 

By Pierre Del Moral, Frederic Patras^ and Sylvain 
rubenthaler 

INRIA Centre Bordeaux Sud Quest, Universite de Nice-Sophia Antipolis 
and Universite de Nice-Sophia Antipolis 

We design exact polynomial expansions of a class of Feynman- 
Kac particle distributions. These expansions are finite and are parametrized 
by coalescent trees and other related combinatorial quantities. The 
accuracy of the expansions at any order is related naturally to the 
number of coalescences of the trees. Our results include an exten- 
sion of the Wick product formula to interacting particle systems. 
They also provide refined nonasymptotic propagation of chaos-type 
properties, as well as sharp Lp-mean error bounds, and laws of large 
numbers for ?7-statistics. 

0. Introduction. The typical phenomenon we are interested in is the 
behavior of particles scattered according to a diffusion process (encoded by 
a Markov transition) and submitted to a potential function. The problem 
is fairly general, since its modelization encodes, besides the scattering of 
particles in an absorbing medium, also filtering problems in signal processing 
or polymerization models, to quote only a few of the applications areas of the 
theory. In these three cases, the potential function represents, respectively, 
the absorption rate, the likehood of the signal values conditional to some 
observations, or the intermolecular attraction or repulsion forces between the 
monomers. For further details and a survey of the applications of Feynman- 
Kac interacting particle models, the interested reader is recommended to 
consult the pair of books [5, 6], and the references therein. 

As it is well known, naive Monte-Carlo methods fail to give a satisfactory 
answer to these questions. This is easily understood with the example of 
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particles in an absorbing medium. In that case, the size of the nonabsorbed 
population will decrease according to the absorption rate of the medium, 
so that, in the end, the empirical distribution of the surviving sample may 
give only a poor approximation of the final repartition of the particles. The 
solution to the decrease of the population size is simple and relies on a mean- 
field approximation. In general, particles will explore the state space as a 
free Markov evolution; during their exploration particles with low potential 
are killed, while the ones with high potential value duplicate. 

During the last two decades, the asymptotic analysis of these Feynman- 
Kac interacting particle models has been developed in various directions, 
including propagation of chaos analysis, Lp-mean error estimates, central 
limit-type theorems, and large deviation principles. The purpose of the 
present work is to improve on these results and develop exact, nonasymp- 
totic, tree-based functional representations of particle block distributions. 
Actually, two probability distributions arise naturally in these models. To 
stick with the diffusion-absorption scheme, they are associated, respectively, 
to the empirical distribution of the nonabsorbed component of the Markov 
chain modeling the trajectories (1.8) and to the corresponding normalized 
empirical occupation measure (1.9). In both cases, the distributions can be 
expanded polynomially with respect to N~'^, where stands for the total 
population size. These Laurent-type expansions are the main results of the 
article. They are described in Theorems 1.5 and 4.12. 

They rely on an original combinatorial, and permutation group analysis on 
a special class of trees and forests that parametrize naturally the trajectories 
of interacting particle systems. Much attention has been paid recently to the 
combinatorics of trees and their applications in physics. A short overview of 
the areas where these phenomena show up is included in Section 3. To the 
best of our knowledge, their introduction in the analysis of Feynman-Kac 
and interacting particle models is new. Since the corresponding geometrical 
and combinatorial study might be useful in other contexts, we have tried to 
isolate its exposition, as far as possible, from its applications to Feynman- 
Kac models. 

The functional expansions allow us to derive as a direct consequence 
Wick formulas for interacting particle systems [Theorem 3.14, (5.1)], re- 
fined, nonasymptotic, propagation of chaos-type properties including strong 
expansions of the particle block distributions with respect to Zolotarev-type 
seminorms (5.2), as well as explicit formulas for Lp-mean error bounds (5.4), 
and laws of large numbers for [/-statistics (5.3), yielding what seems to be 
the first results of this type for mean-field particle models. 

The article is divided into four main parts, devoted, respectively, to the 
precise description of Feynman-Kac particle models and tree-based expan- 
sions (Section 1), to the proof of the expansion formulas (Section 2), to the 
combinatorial and group-theoretical analysis of trajectories in mean-field 
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interacting particle models (Section 3), and to the extension to path-space 
models with applications to propagation of chaos- type properties (Section 4). 
The last section gathers various direct consequences of the main theorems. 

Notation. 

Combinatorial quantities. Let us gather first various notations, for fur- 
ther use in the various parts of the article. 

Let S" be a finite set; we write l^j for the cardinal of S. The notation 
is extended to applications a between finite sets and \a\ denotes then the 
cardinal of the image of a. Precisely, for any pair of integers {l,m) € (N*)^, 
we set [^] = {1, . . . , Z}, and [/][™] the set of mappings a from [m] into [I]. By 
|o|, we denote the cardinality of the set a([m]), and for any 1 <p<l we set 

[;]M = {aG[/]H:|a|=p}. 
When I > m, we also denote by {m,l){= [ifm^) the set of all 

il) 

[l — my. 

one-to-one mappings from [m] into [I], and by Si = {I, I) the symmetric group 
of all permutations of [/]. We also denote by 1^ the identity in S;. Recall that 
the Stirling number of the second kind S{q,p) is the number of partitions of 
[q] into p nonempty subsets, so that 

|[iVf l=S((/,p)(iV)p and N''= ^ S{q,k)iN),. 

l<k<q 

Sequences of maps and associated quantities will be one of the main in- 
gredients of our approach to Feynman— Kac models. We write ^n^q — dcf. 
([gi][9])"+i (resp. IAn,q) for the set of (n + l)-sequences of maps (resp. of 
weakly increasing maps) a = (ap)o<p<n from [q] into itself: 

Notice that we write in bold the symbols for sequences (of maps, integers, 
. . .) such as a. 

For any sequences of integers p = {pk)o<k<n, and 1 = {lk)o<k<n, we write 
p < 1 if and only ii pk < h for all <k <n. We write ||p|| for (pQ + • — \-pn)- 
Assuming now that p < 1, we use the multi-index notation 

n 

(i)p= n (^k)pk, p!= n Pk-^ s{i,p)=Y[s{ik,Pk), 

0<k<n 0<k<n k=0 
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where the s{lk,Pk) are Stirhng numbers of the first kind. Recall that these 
numbers provide the coefficients of the polynomial expansion of the {N)p 
(Stirling formula): 

l<k<p 

see, for example, [3] for further details. The difference (p — 1) and, respec- 
tively, the addition (p + 1) of two sequences is the sequence {pk — lk)o<k<n 
and, respectively, (pk + lk)o<k<n- When no confusions can arise, we write 
N and q, for the constant sequences (A^)o<j<ni and ((/)o<i<n- We also write 
1 and, respectively, 0, for the sequence of unit integers and, respectively, 
null integers. The above definitions are extended to infinite sequence of in- 
tegers p = {pk)k>o S f*^^) with a finite number of strictly positive terms. Any 
function a : N i— > N on the set of integers into itself, is extended to integer 
sequences p S N^, by setting a(p) = {a{pk))k>o- 

In particular, for aG An,q, we write ||a|| the sequence (||ai||)o<i<n- By 
a slight abuse of notation, we write ||a|| for ||(|a|)|| = |ao| + • • • + |a„|. The 
sequence q — |a| is called the coalescence sequence of a. The coalescence 
degree of a is defined by 

n 

coal{a.) := ||q|| — ||a|| = — 

i=0 

The subset of An,q of sequences a such that |a| > q — p is written An,qip)- 

Measures, norms and related notions. Let {E,£) be a measurable space. 
We denote, respectively, by A^(£'), Vi^E) and Bb{E), the set of all finite 
signed measures on {E,£), the convex subset of all probability measures 
and the Banach space of all bounded and measurable functions f on E, 
equipped with the uniform norm ||/|| = sup^g^; |/(x)| . The total variation 
norm is written || • Htv, so that, for any linear operator L on Bh{E), 

||L||tv:= sup \L{f)\. 

/eBf,{E):||/||<l 

A bounded integral operator Q between the measurable spaces {E, £) and 
{F,T), such that, for any / G Bh{F), the functions 

Q{f) -.xeE^ Q{f){x) = Q{x, dy)f{y) G M 

are £^-measurable, and bounded, generates a dual operator fi i— > fj.Q from 
M{E) into M{F), defined by {fiQ){f) := fi{Q{f)). The tensor power Q®^ 
represents the bounded integral operator defined for any / G Bh{F'^) by 

Q®^(/)(xi,...,rE,)= / [Q{xi,dyi)- ■ ■Q{xq,dyq)]f{yi,. . . ,yq). 
JFi 
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For a bounded integral operator Qi from E into F, and an operator Q2 
from F into G, we denote by Q1Q2 the composition operator from E into 
G, defined for any / e Bb{G) by {QMif) := Qi{Q2{f)). 

The notion of differential for sequences of signed measures is also useful. 
Let (0^)7v>i G A^(i?)^ be a uniformly bounded sequence of signed measures 
on a measurable space {E,£), in the sense that supjv>i ||0'^||tv < 00. The 
sequence 0^ is said to converge strongly to some measure @ ^ Ai (E) , as 

I 00, if and only if 

^^fGBbiE) lim e^(/) = e(/). 



Definition 0.1. Let us assume that converges strongly to @. The 
discrete derivative of the sequence {&^)n>i is the sequence of measures 
{de^)N>i defined by 

de^ :=iv[e^-e]. 

We say that is differentiable, if is uniformly bounded, and if it 

strongly converges to some measure dQ £ M{E), as | cxo. 

The discrete derivative of a differentiable sequence can itself be 

differentiable. In this situation, the derivative of the discrete derivative is 
called the second derivative and it is denoted by d^Q = d{dQ), and so on. 

A sequence that is differentiable up to order {k + 1) has the following 
representation: 

0<l<k 

with supjv>i 1 1 9*^^^ 0^1 1 TV < 00, and the convention d^Q = 0, for / = 0. 



1. Feynman Kac semigroups. 



1.1. Definitions. We let {En,£n)n>o be a collection of measurable state 
spaces. In some applications, the series E^ will be constant or there will 
exist canonical isomorphisms between the various state spaces (think to 
the distribution of particles in a fixed absorbing medium), but it may also 
happen that the various state spaces are significantly different (think to 
path estimation and smoothing in signal processing or to polymerization 
sequences of monomers). For a survey of the various applications of the 
theory and examples of families of state spaces En occurring in practice, we 
refer to the pair of books [5, 6], and the references therein. 
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We consider a distribution tjq on £"0, a collection of Markov transitions 
Mn{xn-i,dxn) from En~i into En, and a collection of fn-measurable and 
bounded potential functions G.„ on the state spaces En- We will always 
assume that the potential functions are chosen such that 

(1.1) 0< inf Gn{xn)< sup G„(x„)<oo. 

^r,€En X„&E„ 



We associate to these objects the Feynman-Kac measures defined for any 
function /„ G Bb{En) by the following formulas. 

Definition 1.1. The Feynman-Kac measure 7„ and its normalization 
rjn are defined by 



(1.2) 7n(/n):=IE 



fn{Xn) Yl ^k{Xk) 
0<k<n 



and r7n(/n) :=7n(/n)/7n(l), 



where {Xn)n>o is a Markov chain, taking values in the state spaces {En)n>o, 
with initial distribution r]Q on Eq, and elementary transitions M„ from En^i 
into En- 

By the Markov property and the multiplicative structure of (1.2), we 
check that the flow (r/„)„>o satisfies the following equation: 

(1-3) r/„+i = ^>„+i(r/„), 

where the transformations 'I'n+i '-T^iEn) T^{En+i) are defined for any pair 
iVn,fn+i) G (ViEn) X Sfe(£„+i)) as follows: 

^ f \/f \ r]n{Qn+lifn+l)) 

with Qn+1 ) = Gn{Xn) X Mn+l{Xn,dXn+l). 

The measures ■jn can be expressed in terms of the flow (r/p)o<p<n with the 
formulas 

(1.4) 7„(/„) = 7y„(/„) X 7„(1) with7„(l)= J| r]p{Gp) 

0<p<n 

for any /„ G Bh{En). 

To illustrate their meaning and motivate further the forthcoming devel- 
opments, let us develop with some details an important application of these 
Feynman-Kac flows. Let us consider the distribution of a particle evolving 
in an absorbing medium, with obstacles related to potential functions G„, 
taking values in ]0,1]. In this context, the particle Yn evolves according to 
two separate mechanisms. First, it moves from a site Un-i G En-i, to an- 
other Un £ En according to elementary transitions Mn{yn-i,dyn)- Then, it 
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is absorbed with a probability 1 — GniUn), and placed in an auxiliary ceme- 
tery state Yn = c; otherwise it remains in the same site. If we let T be the 
random absorption time, then it is not difficult to check that 

(1.5) 7n{fn) =nfn(Yn)lT>n] and r]n{fn)=nfn{Yn)\T>n]. 

In the above displayed formulas, we have used the convention that /n(c) = 0, 
when 1^ = c. In time-homogeneous settings (G„ = G and M„ = M), we have 
([5], Proposition 12.4.1) 

7n(l)=IP(r>7^)~e-^'^ 

as n ^ oo. The positive constant A is a measure of the strength and trapping 
effects of the obstacle. By a lemma of Varadhan's, A coincides with the 
logarithmic Lyapunov exponent Aq of the transition operator Q{x,dy) = 
G{x)M{x,dy). Whenever it exists, the corresponding eigenfunction h oi Q 
represents the ground state of the operator Q; see, for example, [5]. 

1.2. Mean-field interacting particle models. A natural mean-field parti- 
cle model associated with the nonlinear Feyman-Kac fiow (1.3) is the E^- 
valued Markov chain sji^^ = {Cn^)i<i<N with elementary transitions defined 
for any F £ Bb{E^) by 

1 ^ 

E(F(e))ieS) = ^n(m(eS))^^(F) with m(eS) = ^E^r- 

J = l 

(1-6) 

In other terms, given the configuration £,^Ji at rank (n — 1), the particle 
system ^i^^ at rank n consists of independent and identically distributed 
random variables with common distribution <I>„(m(^^^^^)). The initial con- 
figuration ^Q^'' consists of independent and identically distributed ran- 
dom variables with distribution r/o. Although the dependency of the model 
on N is strong, due to the mean-field nature of the model, we abbreviate 
{Cn^)i<i<N to ('^^)i<j<Ar when A^ is fixed, except when we want to empha- 
size explicitly the dependency on A^ of the model, for example, as in the 
definition of the measure below. The same observation will apply for 
other random variables, measures or functions in a self-explanatory way. 
Notice that 

^n.(m(ea))(dxj = ^ . fEG„„i(e-_i)M„(e_i,dx„)') 

Z^j=l Lrn-l (,?„_! j \j=i / 

SO that the particle model evolves as a genetic-type model with propor- 
tional selections, and mutation transitions dictated by the pair of potential- 
transition (G„_i,M„). 
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Definition 1.2. The approximation measure 7^ and the normaUzed 
approximation measure 77^ associated with the Feynman-Kac measures 7^ 
and rjn are defined by the empirical occupation measures 

and for any fn£Bb{En), 

(1.7) 7;r(/n):=r/^^(/n)x7^(l) with 7„^(1) := U rj^ (G,). 

0<p<n 

The distribution of particle blocks for any size g < is written : 

(1.8) P^,:=Law(Ci'^,...,er)e^(^^)- 

For a rather complete asymptotic analysis of these measures we again refer 
the reader to [5], and references therein. In particular, for any /„ G Bb{En), 
we have the following almost convergence results: 

lira T]^ (fn) = r]n{fn) and lim 7^(/„) = 7„(/„). 

N—yoo N^oo 

It is well known that 7^ is an unbias approximation measure of 7^, in the 
sense that 

E(7/y(/n))=7n(/n). 

However, it turns out that for any n>l 

E(r?^(/„))=E(/„(e^'^))/r?„(/„). 

This means that the mean-field particle model is not an exact sampling algo- 
rithm of the distributions In practice, it is clearly important to analyze 
the bias of these quantities, and more generally the one of the distribu- 
tions ^nq- From the pure mathematical point of view, the are better 
understood when they are connected with the following measures. 

Definition 1.3. The tensor product occupation measures (resp. the 
restricted tensor product occupation measures) on are given by 



resp. (r?^)®'?:=— - ^ 

^ ae{q,N) 



FEYNMAN-KAC PARTICLE MODELS 9 

The corresponding unnormalized tensor product occupation measures (7^ )®'' 
[resp. restricted unnormalized tensor product occupation measures {'Jn)®'^] 
are defined for any F E I3b{E^) by 

(7^)®nF) := (r,^)^'^(F) x (7^:^(1))'^ and 
(7/r)®^(i^) := {v::f'iF) X (7;^(1))''. 

The corresponding nonnegative measure Q,f^g, indexed by the particle block 
sizes, on the product state spaces E!^, is defined for any F G Bb{E^) by 

(1.9) Qi^{F):=mi:ir'{F))- 

Notice that these various tensor product occupation measures are symmetry- 
invariant by construction. That is, for any F G Bh{E^) and any a G Sg, 

(r?^^r(F) = (r,^r(Foa) and (r?^)®«(F) = o a), 

where a acts by permutation on E^. In particular, we may assume without 
restriction in our forthcoming computations on g-tensor product occupation 
measures, that F is a symmetric function, that is. 

We write from now on, B^f^^{E'^) for the set of all symmetric functions in 
Bb{El). 

Notice also the following symmetry property, essential in view of all the 

forthcoming computations. Since, conditional to ci-ii the are i.i.d., for 
any F G Bb{E^) and any a, 6 G {q, N)'^, we have 

E[F(C(^),...,e:?('^))]=E[F(ef),...,^^))]. 
For instance, we have that 

(1.10) E((r?^)«^(F)) = E{F{en, O) = <,(i^)- 

Let us return to the diffusion-absorption model introduced in the previous 
section. As a consequence of the results presented in [5], the sequence of 
distributions Qn,g converges, as N tends to infinity, to the distribution of q 
nonabsorbed, and independent particles {Y^)i<i<q evolving in the original 
absorbing medium. That is, we have that 

hm Q^^(F) = j^^iF) = E{F{Yl y„^)lri>„ • • • 1t.>„) 

Ai— >oo " — 

for any F G Bb{E^), and where T* stands for the random absorption time 
sequence of the chain {Y^)k>o, with l<i<q. 

This article is mainly concerned with explicit expansions of the deter- 
ministic measures G V{E^), and G Ai{E^), with respect to the 
precision parameter N. 



10 



p. DEL MORAL, F. PATRAS AND S. RUBENTHALER 



1.3. Expansion formulas. Recall that An,q stands for the set of (n + 1)- 
sequences of mappings from [q] into itself. We let A„_g be the nonnegative 
measure-valued functional on An,q defined by 

(1.11) A„,, : a G An,q ^ A^^ := {r^^'^Da^Qf^Da, ■ ■ ■ Q^DaJ e M{El), 

where the operator Db stands for the Markov transition from into itself, 
associated with a mapping b G [q] ^"^1 , and defined by 

A(F)(4,...,x^):=F(xf),...,x^('?)) 

for any F G BbiE^), and (x^ . . . , x^) G E^. Notice that 

DaDb = Dab 

for any pair of mappings (a, b) G ([q']''^')^- For u a linear combination '^i'^i 
of elements of [q] ^''^ , we extend the definition of D by linearity and write 

Du :=^aiDa^. 

Notice that the measures A^ ^ inherit a remarkable invariance property 
from their set-theoretic definition. Namely, let us introduce the natural 
left action of the group S^^^ on An,g defined for all a G An,q and all s = 

(sO,...,Sn+l)GS^+2 by 

s(a) := (soaoS];'\ siois^^, s„anS~| J. 
Then, for any F G Bl^^{E^) we have 
(1-12) Ab^(F) = A^5)(F). 

The identity would not hold if F was not a symmetric function. However, 
as already mentioned, this is not a serious restriction as far as the determi- 
nation of the ^-tensor occupation measures is concerned, since the latter are 
symmetry- invariant . 

Definition 1.4. We write !Fn,q for the set of orbits for the action of 
S^"'"^ on An,q- We write a « b if a and b belong to the same orbit, that is, 
if there exists s G Sg"*"^ such that a = s(b). 

If a G An,q, the orbit of a under the action of s G S^"*"^ is written a. We 
also write Stab (a.) for the stabilizer of a in S^'"^^. According to the class 
formula, the number of elements in a, written #(a) or #(a), is given by 

I btab[a.)\ 
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Explicit formulas for #(a) and for the various quantities associated to the 
action of S^"^^ on An,g such as the number of orbits (and much more) will 
be given later, and, from the combinatorial point of view, form one of the 
cores of the article. 

Notice that, if a ~ b, |a| = |b|, so that the two sequences of maps have 
the same coalescence sequence and coalescence degree. It follows that the 
notions of coalescence degree and coalescence sequence go over to the set 
J^n,q- In particular, notation such as |a| or coal{a) is well defined. In view of 
(1.12), for any f = a in J-'n,q, we set 

K,g = Kg and |f| = |a|. 

The subset of J^n,q associated to An,q{p) is written J^n,q{p)- We are now in 
position to state one of the main results of this article. 

Theorem 1.5. For any 1 < q< N we have the polynomial expansion 

^<k<iq-l){n+l) 

with the collection of signed measures d^Qn,q given by the formula 

r<q-l:||r||=fcfej^„,,(r) '^^''1^1 

for anyF^^B^-^iEl). 

In Section 2, we will prove the exact polynomial expansion (in the pa- 
rameter N~^) of Theorem 1.5. In Section 3, we will develop combinatorial 
concepts in order to compute #(f) for any f G J^n,q (Theorem 3.8). In Sec- 
tion 3, we shall also present several important consequences of Theorem 1.5, 
including explicit descriptions of the first two order terms in the polyno- 
mial expansion, and a new extension of the Wick product formula to forests 
(Theorem 3.14). 

Although more difficult to study since their analysis involves tensor prod- 
uct measures on path-spaces, the distributions have similar expansions 
(Theorem 4.12). The description of these expansions, and of the correspond- 
ing signed measures, is postponed to Section 4. 

2. Particle measures expansions. 

2.1. A preliminary stochastic tensor product formula. The link between 
the two, usual and restricted, tensor products measures, {rj^)®'^ and (??/^)®'', 
relies in the end on a simple observation, that will appear to be fundamental 
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for all our forthcoming computations. Throughout this section, integers N > 
q>l and a measurable state space E are fixed once for all. 
Consider first the surjection 

7T:[qpx{q,N)^[Np, 

{s, a) I — > as := ao s. 

Lemma 2.1. Let be [iV]!"?!; then the cardinal |7r~^({6})| o/7r"^({6}) only 
depends on the cardinal \b\ of the image ofb. It is given by 



|vr-i(W)l = (A^-|fe|),H6|('?) 



\b\- 



Sketch of the proof. Let b£ [A^jt'^l How many possibilities do we 
have to write b in the form b = aos with s G [q]^'^\ a G {q, N)? Set ki,. . . , k^^ G 
[q] such that b{{ki, . . . ,k\^}) = Im{b). Any sequence (without repetitions) 
s{ki), . . . , determines uniquely s (since a is injective and b = ao s) so 
that there are {q)\b\ possibilities for s. Then a(s(fcj)) = b{ki) is fixed for any 
i G {1, . . . , so there remain {N — \b\)g_\^ possibilities for a. □ 

We denote by m(x) the empirical measure associated with an A^-uple 

1 ^ 

1=1 

For any integer q < N,we also consider the empirical measures on defined 
by 

ae[N]M 
^ '1 a&{q,N) 

Although the notation is self-explanatory, notice that, in the sequel, we will 
write simply x'* for {x"-^^\ . . . , x"-''^'>) and 5x« for ^{^'^w ^,,,^x°-{q))- 

Corollary 2.2. We have 

m(x)®^ = m(x)®'?Di^^, 

where 

^1 ~ N't ^ {a\ I 
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Proof. In view of the following identity, that holds for any function 
Fg with g< TV: 

5^a Db{F)=Db{F){x''^^\..., ) = , . . . , x'^''^") ) = 6^at (F) , 

the corollary follows from Lemma 2.1 and the identity 

^ -{N)p{N -p)q_p = l where N>q> p. 



Corollary 2.3. The tensor product measure m(x)'^'' has a Laurent 
expansion: 

m(x)®« = m(x)®''. J2 l^^a^L,, 

0<k<q 

where 

The corollary follows from Corollary 2.2 and from the Stirling formula 

iN)p= s{p,k)NK 

l<k<p 

Assuming from now on that x is generic (i.e., Xj 7^ xj if i 7^ j), we have the 
following corollary. Notice that the assumption that x is generic is here only 
to ensure equalities in Corollary 2.4, if it does not hold, we could replace 
the "=" by "<." 

Corollary 2.4. The following formulas hold for the expansion ofm{x.)'^' 
Af||m(x)®^-m(x)®«||TV q{q-l), 

g 

||m(x)®5Z)gfeiJ|TV = E l^iP^q- k)\S{q,p). 

p=q~k 

The corollary follows from our previous computations, from the definition 
of the Stirling numbers of the second kind and from the following lemma. 

Lemma 2.5. Let Q = J2o<p<q^p'l2ae[N]M (^x", where the Up are arbitrary 
real coefficients; then we have 

iiqiitv= e KmpSiq,p). 

0<p<q 
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Proof. Indeed, we have ||(5||tv < J2o<p<q kpl (^)p5'(g,p) by direct in- 
spection, since ||5x"||tv = 1) for all a G [q]^'^^- Let us write sgn{uk) for the 
sign of Uk, and let us introduce the function cj) G Bb^E'^) defined by 

(Pivi, ■■■,yq)= S9niuk)6^y^ y^^, 

0<k<q 

where A(yi, ...,yq):= \{yi, . . .,yq}\ and (^^(i/i,...,^,,) stands for the Dirac func- 
tion. Then, Q{(j)) =J2o<p<q\''^p\{^)p^iQjP) ™d the lemma follows. □ 

2.2. Laurent expansions. The goal of this section is to prove Theorem 1.5, 
that is, to derive a Laurent expansion of the measures Q^^ G M{E^), Qn,q{F) 
E((7^)®5(F)), F G BbiE^), with respect to the population size parameter 
N. The following proposition is fundamental. Recall from Section 1.3 that 
we write J'n,q for the set of orbits in An.q under the action of the permutation 
group Sg+^. 

Proposition 2.6. For any integers q < N , any time parameter n G N 
and any F G B'^^{E^), we have the Laurent expansion 

n.q ' 

with the mappings A„ ^ introduced in (1.11) properly extended to Tn,q- 

Proof. Combining the definition of the particle model with Corol- 
lary 2.2, we first find that 

= (7/^(1))" xE((,?^f''(F)|^S) 
= (7/^(1))" xE((,?^)0''Z?^^(F)|ea) 
[by (1.6)] = X E(c^„(r?;f_i)«''I)^^(F)|eS). 



Using the fact that, conditionally to and for F G 
and 

iltir'iQ^il)) = (7^-lQn(l))'' = {^:!.,iGn-l)y = (7^(1))", 
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we arrive at 
(2.1) 

Integrating over the past, this yields that 

Using a simple induction, we readily obtain the formulas 

^"^"^'^e^^,, (q)|a| "'"^ ^' 

where An,q is the set of sequences of maps introduced in Section 1.3. In view 
of (1.12), we know that 

a:j,,(f) = aS)(f) 

for any s G S^^^ and the proposition follows. □ 

We are now in position to derive Theorem 1.5. According to Proposi- 
tion 2.6 

<.(^) = ]^ E ^ E #(f)AL(F). 

Using the Stirling formula ([3], page 48), we notice that 

(2.2) (N)p= <vA)N\'\. 

i<i<p 

We also have that 

qUf)= E «(P'i)m7^ E m^l,{F) 

i<i<p<q ^^'^P fe.F„,9,|f|=p 

from which we conclude that 

<.(^)= E E «(P'q-r)]^7^ E #(f)AL(^)- 

r<q-lq-r<p<q V^^P fg^^ _^|f |=p 

The leading term of the development corresponds to the case r = 0. Since 
#(f) = (9!)"^"^ if |f I = *1) this term is equal to 7^*^. Observe then that, for 
r < q — 1, ||r|| < (g — l)(n + 1). The theorem follows. 
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Corollary 2.7. The following formula holds for the expansion ofQ^^^^: 
lim N\\ql^-^T\\^^<{n + l)q{q-l)yi sup 

Remark 2.8. IfO<G<l, then the above supremum is bounded by 1. 
The above corollary is here to show an estimate of the size of the error term 
we derive from Theorem 1.5. 

Proof. The corollary follows from the Laurent expansion of Theo- 
rem 1.5. The term 1/A^ in the expansion is given by sequences r s.t. ||r|| = 1. 
But then, we have two cases to consider, |f | > q — r (which happens if and 
only if |f I = q) and |f | = q — r. 

Then, we have that #(f) = {q\)'^^^ if |f | = q and besides, |{r, ||r|| = 1}| = 
n + 1, and for any such r, s(q, q — r) = — (2) . In conclusion, the case |f | = q 
contributes to (n + V)q{q — l)/2 to the asymptotic evaluation of N\\Q^q — 
7niTV. 

Now, let us consider a sequence r with ||r|| = 1, for example, the one, 
written with = 1, for a given < i < n. In that case, there is a unique 
f with |f I = q — r, which is the set of all sequences of maps a in An,q with 
aj S Sq, j / i and |aj| = q — 1. In particular #(f) = {q\)'^^^{'Q, so that, on 
the whole, the case |f | = contributes to q{q — l)/2 to the evaluation. □ 

To conclude this section, notice that the above functional expansions also 
apply to the dot-tensor product measures 

(7/yr(i^) = (7n^(l)rx(r?^^r(F). 

More precisely, by definition of the particle model we have that 

(2.3) mi^r'{F))=mi:t-ir'QT{F))=Q'i-i,,{QTF). 

3. Combinatorial methods for path enumeration. In the present section, 
we face the problem of computing the cardinals #(f ) involved in our Laurent 
expansions. We also derive various identities relevant for the fine asymptoti- 
cal analysis of mean-field particle interacting models, such as the number of 
elements in J-n.q with a given coalescence degree. To express them, we have 
chosen to use the language of botanical and genealogical trees, both for tech- 
nical reasons (since forests appear to be the most natural parametrization 
of elements in J'n,q)^ and also for the clarity of the exposition. 

Forests, graphs, colored graphs and their combinatorics have appeared 
recently or in various situations in probability and analysis, often in con- 
nection with problems in mathematical physics. A classical example is given 
by the forest formula of Zimmermann, which encodes the renormalization 
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process for the divergent integrals showing up in perturbative quantum field 
theory. The construction is classical, since it goes back to the foundation 
of modern particle physics (see, e.g., [2]), but has been revisited recently 
by Connes and Kreimer, who showed that the algebraic study of Feynman 
graphs (colored by suitable terms appearing in the Lagrangian of the the- 
ory) and of associated trees and forests could give rise to a totally new Hopf 
algebraic approach to the theory [4]. Another example is provided by the 
Butcher group and Runge-Kutta methods in numerical analysis. Here also, 
once again, the fine study of the underlying combinatorial and algebraic 
structures on trees has given rise to a renewal of the subject [1]. Let us 
mention at last the appearance of the combinatorics of colored graphs in 
Gaussian matrix integral models. We refer to [7] and subsequent papers by 
these authors for further references on the subject. 

The approach to Feynman-Kac particle models by means of trees, forests 
and analogous objects which is introduced in the present article is original, 
to our best knowledge, and the methods developed in this setting seem com- 
plementary to the graph and tree-theoretical results that have been obtained 
in other application fields. 

3.1. Jungles, trees and forests. 

3.1.1. Definitions. In this section, we detail the vocabulary of trees, that 
will be of constant use later in the article. The following definitions are 
motivated by the parametrization of tensor product measures expansions 
introduced in the previous sections. 

Let a = (oo, ai, . . . , an) be any sequence of maps : \pk+i] — > \pk], with 
k = 0, . . . ,n. By associating to each element i in [p^] the point {i,k) in 
and for k >0, an oriented arrow from {i,k) to (afc_i(i),/c — 1), the sequence 
can be represented graphically by a graph J(a) in the plane, where the edges 
are allowed to cross. 

More generally, let us consider the set P of graphs in the plane defined 
by a (possibly empty) finite set S of points in (the vertices) and a finite 
subset A of (the arrows or oriented edges) with the following properties: 

(1) The elements of S are of the form (x. A;), € N (notice that, in particular, 
elements in Sk ■= {s S S\3x,s = {x,k)} are naturally left-to-right ordered). 

(2) For each s £ S, there is at most one arrow starting from s. (3) The 
elements of A are of the form ((x, k + 1), (x' , k)), k>0. (4) For any s £ S, 
s = {x,k), there is always a (necessarily unique) sequence of arrows joining 
s to a vertex {x' , 0). 

Two elements {S',A') and {S,A) of P are said to be strongly equivalent if 
and only if there exists a bijection (j) between S and S' with the properties: 
(1) (j) maps bijectively Sk to S"^. (2) (p induces a bijection between A and 
A'. They are said to be LR- equivalent if they are strongly equivalent and 
moreover: (3) (p respects the left-to-right ordering. 
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Definition 3.1. A jungle is an equivalence class in P for the LR- 
equivalence relation. The set of jungles is written J . 

Notice that the process that is described above and maps a sequence of 
maps a to the corresponding jungle [the equivalence class of J (a) in P, still 
written J(a) by a slight abuse of notation] can be inverted. That is, jungles 
are actually in bijection with sequences of maps; however, the introduction 
of jungles is convenient for our purposes, as it will appear below. If j is 
a jungle, we write Seq(X) for the corresponding sequence of maps [so that 
j = J(5'eg(j))]. In the example in Figure 1, we have Seq(X) = (ao)'^i) with cq 
(resp. ai) a map from [7] to [3] (resp. [6] to [7]) and, for example, ai(5) = 7 
and Co (3) = 2. 

Notice that, since the LR-equivalence relation is weaker than the strong 
equivalence relation, the strong equivalence relation goes over from graphs 
in P to jungles. 

Definition 3.2. A forest is an equivalence class of jungles for the strong 
equivalence relation. The set of all forests is written T . A tree is a connected 
forest. The set of all trees is written T. 

Here, connected has the usual topological sense (an oriented graph is 
connected if and only if two arbitrary vertices of the graph can be related 
by a sequence of adjacent unoriented edges). The forest naturally associated 
to a sequence of maps a and to the corresponding jungle J(a) is written 



A tree can be viewed alternatively as an abstract graph — more precisely 
as an abstract finite connected oriented graph without loops such that any 
vertex has at most one outgoing edge. The empty graph is viewed as a tree 
and is called the empty tree. The vertices of a tree without incoming edges 
are called the leaves; the vertices with both incoming edges and an outgoing 
edge are called the internal vertices; the (necessarily unique) vertex without 
outgoing edge is called the root. This terminology, as well as other notions 
introduced below, extend in a straightforward or self-explanatory way from 
trees to forests, jungles, and so on. 

• ^ • ^ • • • • 



F(a). 




\ 




Fig. 1. The graphical representation of a jungle. 
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• • • • • • 




• • • 



Fig. 2. The graphical representation of a planar forest pf = PT1PT2PT1. 

Forests can be viewed as multisets of trees, that is, sets of trees with 
repetitions of the same tree ahowed. An algebraic notation is convenient for 
our purposes, and we write 

f = . . . J'^^k 

for the forest with the tree Tj appearing with multiplicity rrii, i < k. This 
algebraic notation will prove useful, among others, when computing Hilbert 
series for forest enumeration; see Section 3.4. When Tj 7^ Tj for i 7^ j, we say 
that f is written in normal form. 

We say that a tree T is a subtree of T' and we write T CT', if the graph 
of T is a subgraph of the graph of T' , and if the root of T is also the root 
of T'. A subforest f C f ' of a forest f is defined accordingly, as a collection 
of pairwise disjoint subtrees of the trees in f. 

Forests are equivalence classes of jungles. It will be convenient in our 
forthcoming arguments to choose, for a given forest f , a representative in 
J' with nice graphical properties. For that purpose, we introduce the notion 
of planar forest and planar tree. Namely, a jungle j such that Seq{j) is a 
sequence of weakly increasing maps is called a planar forest. A connected 
planar forest is called a planar tree. Graphically, planar forests are jungles 
in which the arrows do not cross. 

It is important to notice that (as illustrated in Figure 2) a planar forest 
pf can be viewed as an ordered sequence of planar trees. Planar forests 
can therefore be represented by noncommutative monomials (or words) on 
the set of planar trees. For example, if PTi and PT2 are two planar trees, 
pf = PT1PT2PT1 is the planar forest obtained by left-to-right concatenation 
of PTi , PT2 and another copy of PTi . 

3.1.2. Functions and operators on jungles, trees and forests. The dis- 
tance between two vertices in a tree is the minimal number of nonoriented 
edges of a path joining them. The height of a vertex is its distance to the 
root. We also say that a vertex with height A: is a vertex at level k in the 
tree. The height ht{T) of a tree T is the maximal distance between a leaf 
and the root. The height ht{i) of a forest f is the maximal height ht{T) of 
the trees T C f in the forest. We write t'fc(f) for the number of vertices in 
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a forest f at level A; > 0. The number of vertices at each level in a forest is 
encoded in the mapping 

v.f eT^v{i) = {vk{i))k>o- 

Notice that for any pair of forests (f , g) G J^, we have, using the algebraic 
notation (the product of forests corresponding to the union of the multisets 
of trees): 

u(fg) =v{f) +v{g). 

We write V for v{J-'), which coincides with the set of integer sequences 
P = {Pk)k>o £ satisfying the following property: 

3/ii(p) e N s.t. inf pfc>0 and sup,,yht(p)Pk = 0. 

k<ht{p) 

For any p G V, we write J-^p for the set of forests f such that v{f) = p. 
Our notation is consistent, since the height ht{{) of a forest f S J-p clearly 
coincides with the height ht{p) of the integer sequence p. When p G V 
is chosen so that po = 1, the set Tp reduces to the set Tp of all trees T 
such that v{T) = p. These notations on trees and forests go over to planar 
forests, planar trees and jungles in a self-explanatory way. The notation is 
also extended to sequences of maps, so that, if p = {po, . . . stands 
for the set of sequences (oq, with Oj a map from p^+i to pi. 

For the unit empty tree Tq = 0, we use the conventions 

v{0) = O={O,O,...), ht{0) = ht{O) = -l, J^o = {0}. 

For any nG(NU{— l}),we denote by V„ C V the subset of sequences p 
such that ht{p) = n. Notice that V_i = {0}. In this notation, the sets of all 
forests and trees with height n are given by the sets 

■■= U -^p and Tn := ]J Tp. 

P&Vn pGVn,PO = l 

The shift operator 

B :p= ipk)k>o G V B{p) = {qi)i>o G V with qi := pi+i 

induces a canonical bijection, still denoted by B, between the set of trees Tp 
with p G V„+i and the set of forests ^B(p) obtained by removing the root 
of the tree: 

B:T€Tp^B{T)€J'B(p)- 

We extend the map B from trees to forests to a map between forests and, 
if f = Tl"' ■ ■ ■ Tf^ is a forest, set 

B{f) = BiTi)"'^ ■ ■ ■ B{Tk)"'K 
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From the point of view of graphs, the operation amounts to removing ah 
the roots and all the edges that have the root as terminal vertex from the 
graph defining f . 

From the graph-theoretic point of view, the number of coalescences c(f)fe 
at each level /c > in a forest f € JF is defined by the mapping 

c:f G^^c(f) ■.= B{v{{)) - |f|. 

The sequence c(f) is called the coalescence sequence of f , following the ter- 
minology introduced for sequences of maps. In the above displayed formula, 
the sequence |f | = (|f|A:)fc>o represents the number of vertices at level k>0 
minus the number of leaves at the same level, that is, the number of vertices 
with an ingoing edge, whereas B{v{{))k stands, by definition, for the number 
of edges between the levels k + 1 and k. Finally observe that for any pair of 
forests (f , g) G we have 

(|fg| =|f| + |g| and B{v{i) + v{g)) = B{v{i)) + B{v{g))) ^ c{fg) = c{f) + c{g). 

We also say that the coalescence order of a vertex in a tree or a forest is its 
number of incoming edges minus 1. The coalescence degree ||c(f ) || of a tree or 
a forest f is the sum of the Cfc(f) or, equivalently, the sum of the coalescence 
orders of its vertices. We say that a tree is trivial if its coalescence degree is 
0. 

3.2. Automorphism groups on jungles. For any given sequence of integers 
p G V„+i, the product permutation group 

Sp := (Spg X Sp-^ X • • • X Sp„_|_j) 

acts naturally on sequences of maps a= (oq, ai, . . . , a„) G ^p, and equiv- 
alently on jungles J(a) G Jp by permutation of the vertices at each level. 
More formally, for any s = (sq, . . . , Sn+i) G Sp this pair of actions is given by 

s(a) := (soaosr\ . . . ,SnOnS~+i) and sJ(a) := J(s(a)). 

An automorphism s G Sp of a given jungle J(a) G Jp is a sequence of per- 
mutations that preserves the jungle, in the sense that sJ(a) = J(a). The set 
of automorphisms of a given jungle J(a) G coincides with the stabilizer 
of a and J(a) with respect to the group action of Sp. 

By definition of forests, two sequences a and b in ^p satisfy F{a) = F{h) 
if and only if J(a) and J(b) differ only by the order of the vertices, that is, 
if a and b belong to the same orbit under the action of Sp. 

Remark 3.3. The set of equivalence classes of jungles in J^p under the 
action of the permutation groups Sp is in bijection with the set of forests 
Tp. In particular, due to the equivalence between the two notions of jungles 
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and sequences of maps, it follows that the set J'n,q of equivalence classes in 
An,q under the natural action of S"^"^^ is canonically in bijection with the 
set of forests with q vertices at levels < i < n + 1 . 

In the next section, we will compute #(f) (the number of jungles asso- 
ciated to a given forest f) by computing the cardinal of the stabilizer of a 
planar forest pf such that -F(pf) =f (such a pf always exists). 

Example 3.4. The jungle j represented in Figure 1 is such that Seq{}) = 
(aO)«i) with 

ai(l) = 2, ai(2) = ai(3) = l, ai(4) = ai(6) = 5, ai(5)=7, 

ao(l)=ao(2) = l, ao(3) = ao(4) = ao(6) = 2, ao(5) = ao(7) = 3. 

For convenience, let us represent a permutation in Sp by the sequence of its 
values on 1, ... ,p, that is, by (ii , . . . , iq) for the permutation s of [q\ such that 
s{j) = ij (take care that this is not the cycle representation of permutations, 
in spite of the notational analogy). We get that the jungle represented in 
Figure 2 is obtained by the action of (123) x (1,234,657) x (312,465) on j. 
The planar forest j' represented in Figure 2 is thus in the same class as j 
(under the action of S(2j,6))- 

A more graphical way of seeing the action of 8(2,7,6) to say that it 
permutes the dots while the edges remain attached to their dots. 

3.3. An inductive method for counting jungles. From the Feynman~Kac 
mean-field approximation point of view, the elements of An,q parametrize 
the trajectories of families of particles. Two trajectories that are equivalent 
under the action of Sq have the same statistical properties; this is the ground 
for the formulas for tensor product occupation measures that have been 
obtained in the previous section. 

In the present section, we face the general problem of computing the 
cardinals #(f), defined as the number of jungles associated to a given forest 
f G J-p, for some integer sequence p G V„_|_i. Let us write f in normal form 
(with the Tj's all distinct) as a commutative monomial of trees: 

f = T™^ • • • r"'' , mi H \-mk=Po. 

Let us also choose arbitrary planar trees PTi representing the T,s, so that 
we can view f as the forest associated to the planar forest 

pf = (PTi)"^! •••(PTfc)™'' 

obtained by left-to-right concatenation of mi copies of PTi, m2 copies of 
PT2 , . . . , rrik copies of PT^ ■ We write PT[ ■ ■ ■ PT^^ for the expansion of the 
noncommutative monomial 

(3.1) pf = (PTi)"'! • • • {PTkT<^ = PT[ ■ ■ ■ PT' 
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as a product of planar trees (without exponents). 
From previous considerations, we have that 

#(f) = # Or6(pf) = #{g G Jp : 3s e Sp s.t. g = s ■ pf}. 

Due to the class formula, we also know that 

#(f) = #o.Mpf) = JMjj = ^, 

where 5ta&(pf) stands for the stabilizer of the jungle pf G i7p- 

The computation of #(f) will be done by induction on (n + 1) = ht{{). Let 
us assume that we know #(g) for any forest of height less than or equal to 
n. Notice that, then, we also know #(t) for any tree t of height {n+ 1), due 
to the canonical bijection B between trees and forests: #(t) = ^{Bt). Prom 
the previous discussion, the problem amounts to computing the cardinals of 
the stabilizers Stab{p{) inductively with respect to the height parameter. 

Lemma 3.5. There is a natural isomorphism 

Stah{{PT^r- . . . {PT,r^)n ({IpJ X S(p,,...,p„^,)) 
~ Stab{PT[) X • • • X Stab{PTl^) 
[with the same notation as in (3.1)]. 

The lemma encodes the observation that an automorphism of a planar 
forest which acts as the identity on the roots decomposes as a product of 
automorphisms of the planar trees in the forest. This follows, for example, 
from the more general fact that an automorphism of a jungle necessarily 
preserves connectedness properties. 

Besides, the first coordinate mapping 

TT : s = (so, . . . , Sn+i) G StabiiPTir^ ■ ■ ■ {PTkT^) 

^ 7r(s) = So G (Smi X • • • X SmJ 

is a surjective map (it even has a natural section, the construction of which 
is omitted), and 

Ker{-K) = ^~Hlpo) ~ Stab{PT[) x • • • x Stab{PT'p^) 

~ StabiPTi)""' X • • • X Stab{PTkTK 

This yields the isomorphism formula 

Stab{{PTir^ ■ ■ ■ {PTkr^)/{Stab{PTir^ X • • • X StabiPTuT'') 

(3.2) 

~ (Srrti X • • • X Smfc)- 

Let us introduce some further notation. 
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Definition 3.6. Let f be a forest written in normal form f = T{"^ • • • T™*" 
In this situation, we say that the unordered /c-uplet (mi, . . . , m^) is the sym- 
metry multiset of the tree T = B~^{i) and write 

s(r) = s{B-\Tr ■ • -Tf *)) = (mi, . . . ,mfc). 

The symmetry multiset of the forest f = T™^ • • • T^*" is the disjoint union of 
the symmetry multisets of its trees 

s(rr • . .T^'-) = s(ri), . . . ,s(ri) , . . . xn), . ..,s{Tk) ] ■ 

\ rni terms terms / 



For example, the symmetry multiset of the first tree of the forest f dis- 
played in Figure 2 is (1,1); the symmetry multiset of B~^{{) is (2,1); and 
the symmetry multiset of f is (1, 1, 3, 1, 1). 

The definition of symmetry multisets extends naturally to planar forests 
and planar trees. In particular, if pf is a planar forest with f as its underlying 
forest, we set s(pf) := s(f). 

We deduce from (3.2) the following recursive formula for the computation 
of #(f). 

Proposition 3.7. We have 

k 

I StabiiPTir^ ■ ■ ■ {Pnr^)\ = n(m.!| Stah{B{PT{))r^) 

i=l 

= s{B~\jyi))\X{{\Stah{B{PT,))r). 

1=1 

Theorem 3.8. The number of jungles (or, equivalently, of Feynman- 
Kac mean-field type trajectories with the same statistics) in i € J^p, with 
p G V„+i, is given by 

p! 

where we use the usual multi-index notation to define s(i?*(f))!. 
Proof. It clearly suffices to prove that 

n 

\Stabip{)\= n s(i3*(pf))! 

i=-i 
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for any planar forest pf of height (n + 1). We check this assertion by induc- 
tion on the height parameter. First, we observe that a planar forest pf of 
height 1 can be represented as a noncommutative monomial 

pf = (PTi)™! •••(PTfcr'^ 

with different planar trees PTi of height 1, and some sequence of integers 
TTij. In that case, the planar forests B[PTi) reduce to sequences of trees with 
null height. This yields that 

\Stah[B{PTi))\=s{PT{)\. 

By Proposition 3.7, we conclude that 

k 

I 5to6(pf)| = \{{m,\{s{PT,)\r^) = s(S-i(pf))!s(pf)!. 

This ends the proof of the formula at rank 1. 

Suppose now that the assertion is satisfied for any planar forest with 
height at most n. By Proposition 3.7, for any planar forest pf = [PTi)"^^ ■ ■ ■ 
{PTk)^^ , with height (n + 1), and written in terms of distinct planar trees 
PTi, we have 

k 

I 5te6(pf)| = s(p-i(pf))! \{ I Stab{B{PTi))r\ 

1=1 

Since the planar forests B{PTi) have height at most n, the induction hy- 
pothesis implies that 

n— 1 n 

I Stab{B{PTi))\ = n ^{B\B{PTi)))\ = [] s{B\PT,))\. 

j=-l j=0 

Recalling that -B^ (pf) = Bi{PTi)'^^ ■ ■ ■ B^ {PTk)"'^^ , we also find that 
s{B^{p{)) = (s(i?^(Pri)), . . . ,s(i?^(PTi)), . . . , s{B^{PTk)), . . .,s{B^iPn))). 

" V ' V ' 

mi terms mj. terms 

We conclude that 

n 

I Stabipi)\ = siB-\p{))l n s(i?^(pf))!. 
The theorem follows. □ 
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3.4. Hilbert series method for forest enumeration. In the present sub- 
section, we face the problem of computing other cardinals relevant to the 
analysis of Feynman-Kac mean-field particle models, according to our The- 
orem 1.5. For example, we want to be able to compute the number of forests 
in any JTp [the set of forests f such that u(f) = p], or the number of forests 
in J^n.q (as defined in Definition 1.4) with a given coalescence degree or a 
prescribed coalescence sequence. Recall that the notions of coalescence se- 
quence and degree, as defined in the Introduction for sequences of maps in 
A riM, go over to arbitrary forests and trees. 

Recall first some general properties of free monoids. Let us write {S) for 
the free commutative monoid generated by a set S: 

(S) = {/f 1 • • • : A; > and VI < i < fc, > and fi G 5} 

with the convention /{"^ • • •/^'■" = 1, for k = 0. We then have the Hilbert 
series (i.e., formal series) expansions 

(3.3) r^=E-" nT^= E ^• 

Now, if ^ is a multiplicative map from {S) to an algebra A, we have, whenever 
the expression on the right-hand side makes sense (this will be obviously the 
case in the examples that we will consider): 

Since a forest is, in the algebraic representation, nothing but a commuta- 
tive monomial of trees, this result applies to forest enumeration. Moreover, 
as we shall see, suitable refinements of the previous identity lead naturally to 
a formal series enumeration approach to the various quantities meaningful 
in the study of Feynman-Kac particle models. 

Definition 3.9. The (multidegree) Hilbert series of forests, EI^(x), is 
the Hilbert series associated to the partition of the set of forests of height 
less than or equal to n into subsets according to the number of vertices at 
each level, that is, 

H^(x):= E l-^p|xP=M^-i(x)+ ^ |.Fp|xP, 

ht{p)<n peVn 

where we write x^ shorthand for 

We let dp = • • • , and we consider the mapping B from V into 

itself defined for any /c > by 

S-i : p e Vfc ^ B-\p) = (1, p) = (l,po,Pi, ■ ■ ■) e Vfc+1. 
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Since forests f of height are characterized by the number of roots in f , we 
clearly have the formula 

(Vp G Vol^pl = 1) ^ M%{xo) = = ^ xg. 

1 — Xn ^ 



Proposition 3.10. 

I^pI 



Vn>l ]HlKx)= n ( i_,b-Mp) ) 



ape"-i(o) 



pev„ 

/n addition, the generating function associated with the number of forests 
with multidegree p G V„, and with trees with height n (elements ofTn), 
is given by 

/ ^ xapH^-i(o) 
Vn>l 5: KT„)n.Fp|xP= n L_,b-Hp) 



The proposition follows directly from the formula for Hilbert series of free 
commutative monoids and from the following three observations. First, the 
map that sends a forest in !Fp to is a multiplicative map from the set 
of forests (viewed as a free commutative monoid) to the algebra of formal 
power series over the variables Xj. Second, as we have already observed, the 
set of forests identifies naturally with the free commutative monoid over 
the set of trees, another formulation of the fact that a forest is a multiset 
of trees. Third, the maps and B define bijections between the set of 
forests of height n — 1 and trees of height n. 

Notice that these formulas make the Hilbert series computable at any 
finite vertex order, and at any height using any formal computation soft- 
ware. The first two orders can be handily computed. We already know the 
generating series for n = 0. By Definition 3.9 we find 

]HI^(x) = 5: + \^pK^i'- 
Po>o peVi 

This readily yields that for any p € Vi , we have 

l-^pl = # S {K)n>o G : Ykn=Po and ^ nkn =Pi>- 

[ n>0 n>0 J 
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More generally, recall that for any m > 1 , we have the following formal series 
expansion 

1 \"' _^{m-l + k)lf^ 



so that 



irn+l 



(x)= n E 



u 



l-^pl - l + ^p)! „fc„B-i(p) 



W - 11 (IT \- l)\k ' 



^ ( Y[ (l-^pl-i + Mp))! V^y:..„..„^-(P)B-^(P)_ 

fe:{p:fti(p)<n}^N \fei(p)<n d-^P I ~ ^) '^(P) ' / 

This implies that for any p € V^+i, we have 

fc:{q:M(q)<n}^N \ht{q)<n ^^'^'i^ ^J-^^W' / 

E.U<,)<n'=('l)"^-'('l)=P 

To compute or estimate the Laurent expansions of the distributions Q^g, 
we are actually interested in a more precise Hilbert series, namely the one 
taking into account, besides the multidegrees of forests, their coalescence 
numbers. 

Definition 3.11. We denote by -Fp[q] := J"p n c~^(q) the set of forests 
in J^p with a prescribed coalescence sequence q (we also use the convention 

We let C^(x,y) be the (multidegree) Hilbert series of forests associated 
to the partition of the set of forests of height less than or equal to n with 
prescribed coalescence sequences, that is, 

c^(x,y):= J2 E l-^p[q]|xPy'^- 

ht{p)<n^eC{p) 

In the above display, C(p) stands for the set of coalescence multidegrees: 
Vp G V C(p) := c(.Fp) = {q G : {B{p) - p)+ < q < {B{p) - 1)+} 
with (-B(p) - 1)+ = {{pk - l) + )fc>i, for any p = {pk)k>o e V. 

Proposition 3.12. The multidegree Hilbert series of coales cent forests 
C^(x,y) satisfies the recursive formula 

C^(x,y)=Cr^(x,y) Y[ Y[ ( . / 

peV„_iqGC(B-i(p)) ^ ' 
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Notice first that JFp = UqGC(p) -^pMi ^^^d ^t{q) < {ht{p) — 1), for any 
q G C (p) . Also observe that 

^b-Hp)=T^i,p) =^ C(i?-Hp))=c(r(i,p)) = {po-l}xC(p). 

Notice at last that the map from forests to formal power series that maps a 
forest in ^p[q] to x^y^^ induces a multiplicative map from forests to formal 
power series. The proposition follows then once again from the formula for 
Hilbert series of free commutative monoids and the correspondence between 
trees and forests. 

Using the same lines of arguments as before, this proposition can be used 
to derive a recursive formula for the explicit combinatorial calculation of 
the number of forests with prescribed heights and coalescence multi-indices. 
The first two orders are given by 

Ci(x,y) = and Cl(x,y) = fT — . 

Notice that the Hilbert series technique can be developed to any order of 
refinement. For example, it could be used to take into account, besides the 
number of vertices or of coalescences at each level, the cardinals #(f), so that 
the coefficients of the Laurent expansion of the measures could be read, 
in the end, entirely on the corresponding Hilbert series. We leave the task of 
expressing the recursive formula to the interested reader, and simply point 
out that the technique allows an easy, systematic, recursive computation of 
the coefficients of the expansion of the measures Qn,q order, in both 

n and q. The observation can be useful, especially in view of the systematic 
development of numerical schemes and numerical approximations based on 
Feynman-Kac particle models. 

3.5. Some forests expansions. In the present section, and the forthcom- 
ing one, we take advantage of the language of trees and of the results ob- 
tained on their statistics to compute the first orders of the Laurent functional 
representation of Qn,q and, respectively, to derive a natural generalization 
to Feynman-Kac particle models of the classical Wick product formula. 

Let us fix n and (/ > 4, so that the notation q denotes, once again, the 
constant sequence of length n associated to q. As we have already pointed 
out, there is only one forest in J^^ without coalescence, which is the product 
of q trivial trees of height n. There is only one forest in J^^ with only one 
coalescence, at level i, that will be written fi^j. Its only nontrivial tree is the 
tree with one coalescence at level i and two leaves at level n. There are two 
forests with coalescence degree 2 and two coalescences at level i, written 
£2 j and fl j. The notation j denotes the forest with only one nontrivial 
tree with coalescence degree 2, a vertex with coalescence order 2 at level 
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Fig. 3. Examples of planar tree representatives pfJS ^ of the f^ ,, 's for q — 4, n = 2. 



i and its three leaves at level n. The notation denotes the forest with 



two nontrivial trees with coalescence degree 1 and the coalescence at level i. 
There are four forests with coalescence degree 2 and two coalescences at 
levels i <j, written f2ij; A; = 1,...,4. The forest fgjj has one nontrivial 
tree with coalescences at levels i and j and its tree leaves at level n. The 
forest fl j has one nontrivial tree with coalescences at levels i and j and 
its three leaves at the levels j,n,n. The forest ff j j has two nontrivial trees 
with one coalescence at level i, resp. j and their two leaves at level n. The 
forest f|jj has two nontrivial trees with one coalescence at level i, resp. j 
and their two respective leaves at level j,n, resp. n,n. See Figure 3. 

Expanding the formulas for 9*Qg,nj using the formulas obtained in Theo- 
rem 3.8 for the cardinals #(f ) and using that s{q, q — 2) = (3) (see, e.g., 
[3], on page 63), we get the following result. 



Corollary 3.13. The first three order terms in the polynomial func- 
tional representation of Qn,q ^'^^ given by the following formulas: 

d'Qn,,= '-^^ E (<'|-7r), 
0<fc<n 

4 y 

^ ^ 0<k<l<n 

H 2^ ( A„;g' +(g-2)(A„;g +A„,g') 

0<k<l<n ^ 
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, (g-2)(g-3) ^f3,,^ 

3.6. A Wick product formula on forests. Let now Bq^^{E^) GBb{E^) be 
the set of symmetric functions F on such that 



{Dl^-l 'S>'yn){F){xi,...,Xq^i) = J F{XI, . . . ,Xg^i,Xq)-fn{dXq) =0. 

Notice that B^'^{E'i[) contains the set of functions F = (P)synn with P £ 
Poly(£'^), where Poly(£^^) C Bb{E^) stands for the subset of polynomial 
functions of the form 

P = Y^ c{a)r with = (^(1) ® • • • ® 

In the above display, / is a finite subset of N^''\ c G M^, and the elementary 
functions E Bb{En) are chosen such that 7n(/^) = 0. For instance, we can 
take 

r=(/-^?n(5^)) with/ES,(S„). 
Notice first that, by Theorem 1.5, for any integer q and any F G Bq^^{Eq), 



and, if q is even. 



29/2(g/2)! 



where f is the forest in .Fo,q containing | copies of the tree of unit height, with 
two vertices at level 1, and | copies of the tree with the root as unique vertex 
(the reader may check these identities or refer to the more general arguments 
given below). For symmetric tensor product functions F = (/^ ® • • • f'^)sym, 
associated with a collection of functions /* G Bb{Eo), such that 77o(/*) = 
7o(/*) = 0, for any 1 < f < q^, we readily find that 

(3-5) AS,,(^) = ^^^%^Ef n ^o(/7^)) 

as soon as q is even. In the above displayed formula, Xq ranges over all par- 
titions of [g] into pairs. In a more probabilistic language, the above formula 
can be interpreted as the gth-order central moment 

(3.6) 

^'^M^^Y.l n E(W^o(r)T^o(/^))) 
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of a Gaussian field Wq on the Banach space of functions BbiEo), such that for 
any pair of functions (93, V') G Bb{Eo), E{Wo{ip)) = and E{Wo{^)Wo{ip)) = 
r]Q{(p'ip). In the above equation (3.6), the second equahty is the classical Wick 
formula. 

We choose to call (3.5) and the more general identities that will appear 
below "Wick formulas" because of this interpretation. More precisely,in the 
present section we will show that the Wick formula (3.5) generalizes to 
forests in Tn,q of arbitrary height. Let us start by listing various straight- 
forward properties of trees and forests. A tree T with coalescence degree d 
has exactly {d+ 1) leaves. A forest with coalescence degree d has at most d 
nontrivial trees, and the equality holds if and only if all its nontrivial trees 
have coalescence degree 1. In particular, if a forest in J-n,q has coalescence 
degree d, it has at most {2d) leaves belonging to nontrivial trees so that, if 
(i < |, there is at least one vertex at level n + 1 belonging to a trivial tree 
(i.e., a tree with coalescence number 0). 

The same reasoning shows that, when d = |, a forest in J-'n.q with coales- 
cence degree d, and coalescence sequence r, does not contain a trivial tree 
of height (n + 1) if and only if it is the forest fr := Tq°Uq° ■ ■ ■ T^"Un" , where 
we write for the unique tree of coalescence degree 1 with a coalescence 
at level k and its two leaves at level (n + 1), and where we write Uk for the 
trivial tree of height k. 

We conclude this series of remarks by noting that, if f G J^n.q can be 
written as the product (or disjoint union) of a forest g in Tn,q-i with Un+i, 
the trivial tree of height (n + 1), then, for any F G B^^'^{E^) we have, by 
definition of the measures 

which is equal to if F G B''^'°{El). 

We are now in position to derive the forest Wick formula using Theo- 
rem 1.5. 



Theorem 3.14. For any even integer q < N and any symmetric func- 
tion F G B^"^{E'^), we have 

yk<q/2 d''QnqiF)=0 and 
(3.7) ' , 



r<q,||r|| = | 



For odd integers q < N , the partial measure-valued derivatives are the 
null measure on Bq^^{E^), up to any order k < [q/2\. 
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We close this section with a Gaussian field interpretation of the Wick 
formula (3.7). We further assume that q is an even integer. We consider a 
collection of independent Gaussian fields {Wk)Q<k<n on the Banach spaces 
{Bb{Ek))o<k<n, with for any {ipk,ipk) G Bh{Ek)'^, and < A: < n 

]E(VFfc(^fc)) = and ¥.{Wk{^k)Wk{ilJk)) = lk{^k^k)- 

We also introduce the centered Gaussian field Vn on Bb{En) defined for any 
If G Bb{En) by the following formula: 



^nM= E yik{l)Wk{Qk,n{^))- 
0<k<n 

Let {(pi)i<i<q G BbiEnY be a collection of functions such that Jn{fi) = 0, 
for any I <i <q. For the tensor product function F = J2aeSq {^cr{i) C5 • • • ^ 
93o-(g)), one can check that 

At(^)=^E n T.Wiiy n ik{QkAvi)QkAvj))\- 

In the above displayed formula, the first sum is over the set I of all partitions 
/ = {Ik)o<k<n of [q] into (n + 1) blocks with cardinality \Ik\ = (2rfc), and the 
second sum ranges over the set Ik of all partitions Jk of the sets Ik into 
pairs, with < A; < n. By definition of the Gaussian fields {Wk)o<k<n, and 
due to the classical Wick formula, we find that 



2<?/2r' 

l£XO<k<n \ielk 



U I[^Jlkil)WkiQk,n{^^)) 

\IeIO<k<niGlk ' 



from which we arrive at 



E lS^^^nAF.^)=A E E n n\/^^'^(QM(¥'.))). 

r<q,||r||=g/2 ' Vr:||2r||=i} / 0<k<ni€lk / 

Recalling that all Gaussian fields {Wk)o<k<n are independent and centered, 
we get 

5^/'Qn,,(i^n)=E( n ( E ^/liJ^)WkiQkAV^))]] =n R ■ 
\l<i<q\0<k<n // Vl<i<g / 

Written in a more synthetic way, we have proved the following formula: 
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This result can alternatively be derived combining the Lg-mean error es- 
timates presented in [5], Theorem 7.4.2, with the multidimensional central 
limit theorems presented in [5], Proposition 9.4.1. More precisely, the q- 
dimensional particle random fields 

converge in law to {Vn{<fn))i<i<q- By the continuous mapping theorem, com- 
bined with simple uniform integrability arguments, one checks that 

4. Extension to path-space models. In the present section, we extend 
our previous analysis to the statistical study of path-spaces. Our goal is 
now to prove Theorem 4.12 which states an expansion for measures 
(analog to Theorem 1.5 for measures Qn,g)- It will appear in the proof of 
Theorem 4.12 that we need expansions of tensor products of measures (of 
the type 7^*^^ ® • • • (8)7^'^"). To this purpose, we will introduce and study (in 
Sections 4.1-4.3) colored trees and forests. We will also define Feynman-Kac 
semigroup on tensor products of measures (4.4) and finally derive expansions 
of tensor products of measures (in Section 4.5). Unfortunately, Theorem 4.12 
can only be written down after all these technicalities. 

Let us make two general remarks about colored forests. First, the (two) 
colors are used to distinguish between vertices that are authorized or not 
to have an incoming edge. Second, the statistical meaning of colored forests 
is similar to the one of forests as far as the approximation measures 7^ 
are concerned. Namely, the introduction of colors on vertices is a natural 
consequence of the replacement of 7^ by the path-space measures {'^0)^'^° 
■ ■ ■ ® (.In)^'^" ■ However, a precise understanding of the meaning of colors can 
hardly be obtained without entering the details of the constructions, and we 
refer therefore to Section 4.5 and Theorem 4.2 for deeper insights. 

4.1. Colored trees, forests and jungles. 

Definition 4.1. A colored tree is a tree in T with colored vertices, 
with two distinguished colors, say black and white. Only black vertices may 
have an ingoing edge. That is, equivalently, all white vertices are leaves — the 
converse being not true in general. A colored forest is a multiset of colored 
trees. The sets of colored trees and colored forests are denoted, respectively, 
by T and J-'. 

As usual, a colored forest can also be viewed as a commutative monomial 
over the set of colored trees. In this interpretation, the generating series 
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techniques that we have developed to deal with the enumeration of forests 
will apply to colored forests. The computation of the corresponding series is 
left to the interested reader. 

Most of the notions associated to trees and forests go over in a straight- 
forward way to colored forests and colored trees. As a general rule, we will 
write a line over symbols associated to colored trees and colored forests. For 
instance, we write Vk(T) = {vk{T),v'f^(T)) for the number Vk{T), respectively 
v'ij.{T), of white, respectively black, vertices in the colored tree T, at each 
level k>0. _ 

We also let V be the set of all sequences of pairs of integers p = {Pk)k>o £ 
(N^)^, with pf. = {pkiP'k) for every > 0, satisfying the following property: 
3M(p) G N s.t. 

inf pk + Pfc > 0> P'k > and sup Pfc + Pfc = 0. 

k<ht{p) k<ht{p) k>kt{p) 

For any n G N, we denote by V„ C V the subset of sequences p such that the 
height of p, ht{p), is equal to n. Finally, we let Tp be the set of colored trees 
T with Vk{T) =Pk white vertices, and v'^{T) =p'^ black vertices, at each level 
/c > 0. Since T is a colored tree, this implies that po = and = 1, except 
if /ii(p) = 0. In that case, po m^-Y also be equal to 1 (and then = 0). 

Definition 4.2. For p G V^+i, let us cah, by analogy with the un- 
colored case, colored jungles (resp. colored planar forests) of type p any 
sequence of maps (resp. of weakly increasing maps) (oq, • • • , a„), where 



The set of colored jungles is written J . 

The permutation group analysis of jungles and forests can be extended 
to the colored case. For that purpose, we write 



the group is acting on [p] by (a, /3)(i) := a{i) and on \p'\ by {a,f3){j) := I3{j). 

We can define colored planar forests in a graphic way by saying they are 
colored jungles such that, at any level, edges starting from a fixed color 



ak = iak,ai)G[piyP^+^^ x [p'^] 




o 




Fig. 4. A colored jungle which is a colored planar forest. 
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(black or white) do not cross each other. It is the case in Figure 4, although 
it should be noticed that some vertices are crossing each other (so that the 
notion of planarity for colored forests is not the intuitive one). 

Once again, the notation we use on colored forests extends to colored 
jungles in a self-explanatory way. For instance, J^, p G V„_|_i, stands for 
the set of colored jungles j with Vk{i) = Pk white vertices, and t'fc(j) = p'k 
black vertices, at each level < A; < (n + 1). For any p G V„_|_i, we also have 
that Jp ~:4p, where :4p := [p'q]^p^^ x x • • • x [p'JP"+'^ x bn]'^"+i^- 

4.2. Automorphism groups on colored jungles. As for the expansion of 
measures in Section 2, we will need a definition of equivalence classes of 
colored jungles, a parametrization of the classes and a computation of their 
cardinals. It appears that most of the constructions in Sections 2 and 3 can 
be extended without further work to colored jungles. Since all constructions 
can be mimicked, we only outline the main ideas of the generalization from 
jungles to colored jungles. 

We let p G Vn+i be a sequence of (n + 1) integer pairs = {pk,p'k), with 
< A: < (n + 1). We associate with p the product permutation group 

Sp = Sp^ X Spj X • • • X Sp^^^ . 

This group acts naturally on sequences of maps a = (ao,ai, ...,«„) G -4p, 
where a, G x [p^j'^'^+i^, and on jungles J(a) in by permutation of 

the colored vertices at each level. More formally, for any s = (sq, . . . , Sn+i) G 
Sp this pair of actions is given by 

s(a) := {-soaQsl^, . . . ,s„a„s^:{:J and sJ(a) := J(s(a)). 

Two colored jungles in the same orbit under the action of Sp have the same 
underlying colored forest. Conversely, if two colored jungles in have the 
same underlying colored forests, they differ only by a permutation of the 
vertices of their colored graphs that preserves the colors of the vertices, and 
therefore are in the same orbit under the action of Sp. In other terms, we 
have the following lemma. 

Lemma 4.3. Equivalence classes of colored jungles in J-p under the ac- 
tion of the permutation groups Sp are in bijection with colored forests in 
^p. 

As for usual trees and forests, we write B for the map from colored trees 
to colored forests defined by removing the root of a colored tree, and we 
write also, as in Section 3.2, B for the induced map from the set of colored 
forests into itself. 
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4.3. An inductive method for counting colored jungles. Let us conclude 
this section by enumerating the number #(f ) of colored jungles associated to 
a given colored forest f G JFp, with p G Vn+i - The process is as in Section 3.2, 
and the result follows ultimately from the class formula 



#(f) 



Stab{i)\ 

In the above displayed formula, p! stands for the multi-index factorial p! = 
rifcioPfc'Pfc'! s-i^cl we have written abusively \Stab{f)\, for the cardinal of 
the stabilizer in Sp of any representative pf of f , where f is viewed as an 
equivalence class of colored jungles. 

Let us assume that f can be written, as a monomial over the set of colored 
trees, as 

f=n^' 



=1 



where the Tj's are pairwise distinct and mi > • • • > m^. We write PTi for 
a set of representatives of the Tj viewed as equivalence classes of colored 
jungles (beware that, here, the prime exponent has no color meaning). 

As in Section 3.3, we shall write s(i?~^(f)) for the unordered fc-uplet 
(m-i, . . . ,m/c), called the symmetry multiset of the colored tree B~^(i), and 
extend the notation to colored forests so that s(f ) is the disjoint union of the 
s(Tj). The proof in Section 3.3 goes then over without changes, except for 
the replacement of forests, planar forests and planar trees by their colored 
analogues. For instance, we have the recursion formula 

k 

\ Stab(^)\ = Y[{mi\\ Stab B(PT^)r^). 
1=1 

Specifying the recursion, we deduce the following theorem. 

Theorem 4.4. The number of colored jungles with underlying colored 
forest f G J-p, with p G Vn+i, is given by 



nr=_is(i?Hf))!' 

where we use the usual multi-index notation to define s(i3*(f))!. 

4.4. Feynman-Kac semigroups. In this section, the parameter n > rep- 
resents a fixed time horizon. For any sequence of integers q = (qq, . . . ,qn) G 
N""*"^, and any — 1 < m < n, we set 

(4.1) q',,:= J2 9k. 

m<k<n 
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Notice that, for m<n, c(m — Qm.+i + Qm+i- We associate with q the unnor- 
mahzed Feynman-Kac measures 

7^ = • • • T®*?" G (E^) where B^ = E^q°x---x E^^ 

is equipped with the tensor product sigma field. Points in are indexed 
as follows: 

((xg, ■ • ■ , Xq*)-, ■ ■ ■ , (Xn, • ■ • ) 3;^")). 

Definition 4.5. We let (7p')o<p<n, and (Q^)o<p<n, be the collection of 
measures and integral operators defined by 

with the operators qI^'' I'^p-i) ^^^^ e'^_^^'^^^^ into {E'^^Si x Ep'^^) defined 
by the tensor product formula 

(recall that the 7p's and Qp's have been defined in Section 1). 

Notice that, for any < p < n, 7^" is a positive measure on the product 
space 

E^ := X • • • X E'/_\' X E'p'^"'^ 
and Qpf is a positive integral operator from E^^^^ into E^. 

Lemma 4.6. For any sequence of integers q G N"~^^, and any <p<n, 
we have 

In ip ^p.n 

with the semigroup {Qpj^^p^)i<pi<p2<n defined by 

Qpi,P2 ~ Qpi+i(Qpi+i,P2) (~ Qpi+i ■ ■ ■ Qp2"iQp2) 
and the convention Q^f^^p^ = W, the identity operator, for pi =P2- 

Proof. We start with observing that 

Therefore, we find that 7^ = 7^_]^Q5^. Using a simple induction, the lemma 
follows. □ 
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4.5. Unnormalized particle measures. In this section, we derive a func- 
tional representation and a Laurent expansion of particle tensor product 
measures on the path-space similar to the ones obtained in Theorem 1.5. 

In the further development of this section, the time horizon n > is a 
fixed parameter, and we let q = (go, . • . ,qn) G N""*"^ be a given sequence of 
{n+ 1) integers. Before we turn to the study of the particle tensor product 
measures 

and the associated path-space measures Q^q G A^(E^) defined by 

let us introduce some useful definitions. We associate with q the pair of 
sequences q and q' defined by 

q := {Qk,qk)-i<k<n and q' := {q'k)^i<k<n 

with the convention q^i = 0, and the integer sequence q' introduced in (4.1). 
Notice that q'^ = 0, ||q|| = qo + q'o, ht{q) = {n + 1), and 

|q'l= E 9fc= E iik + (i'k)= E + 

-l<fe<n. 0<fc<n 0<fc<n 

Notice that [g^.j'^^-i^ = bfc„i]'^''l x [g^^j'^fcl. We will identify the two 
sets of maps throughout all the sequel. 

Definition 4.7. For G Wk^i]^'^''^ x [q'k~i\^'^'''K we define the Markov 
transitions 

For u a linear combination X^ie/CKjOfc^ °^ ^i* ^ ^ fefc-il^'^'^^ we ex- 

tend the definition of D by linearity and write 



D^,. = E«^DL( 



We let A„^q be the nonnegative measure- valued functional on defined by 
An,q:a' = (ao,ai, . . . ,a„) 

^ A^;q := ryo^l'^lD^-^Q^D?- • • • Q^D;^-^ G M{K)- 

Definition 4.8. We write ""(E^) for the elements F of Sb(E^) that 
are symmetric in the first q^ variables, the next qi variables, . . . , and the 
last qn variables. That is, we have that 

V(ao,...,cT„) G (Sgo X ••• X SgJ F= (Z?,„ •••®L',JF. 
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We consider the subset i3o^°'(E^) C "(E^) of all functions F such that 

for any <p<n. 

Notice that B^^(E^) contains the tensor product 

{B'J'^iE'o') ® • • • ® So"iE?r)) = {Fo ® • • • ® i^n, : VO < p < nFp G B'J'^iE^/)} 

of the sets Bq^"^{Ep') introduced in the beginning of Section 3.6, with < p < 
n. For the same reasons as in (1.12), we also have the symmetry invariance 
property 

VF G B^iE^) , Vs G Sq, Vj G AJ;q(F) = A^g (^f) , 

so that q is well defined for f a colored forest. By construction, we also 
have that 

as soon as the colored forest f contains at least one trivial colored tree with 
a white leaf. 

Recall that a tree or a colored tree is said to be trivial if its coalescence 
sequence is the null sequence of integers. Notice also that for any 1 < q', the 
set of colored forests with exactly lk{& [0,(?fe_i[) coalescent edges at level k 
is given by 

:={fG:Fq:|f|=q'-l}. 

Thus, for any r < q', the set of colored forests .Fq(r) with less than r/j 
coalescent edges at level k is given by 

7=q(r):=U7=q[l]. 

l<r 

The coalescence degree of a colored forest f G .^q[l] is the sum |q' — 1| of the 
coalescence orders of its vertices. 

Notice that a colored forest with coalescence degree d has at most d non- 
trivial colored trees. In addition, a colored forest with a coalescence degree 
d has at most {2d) leaves belonging to nontrivial trees. Since a colored forest 
in jTq has exactly ||q|| white leaves, if ||q|| > {2d), then it contains at least 
one trivial colored tree with a white leaf. 

Next, we discuss the situation where ||q|| is an even integer, and we char- 
acterize the subset of forests (in ^q), with a coalescence degree d= ||q||/2, 
without any trivial colored tree with a white leaf. This characterization fol- 
lows the same lines of arguments as the ones presented on page 808. For any 
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integers 0<k<l<m<n, we let Tk,i,m be the unique colored tree, with 
a single coalescence at level k, and white leaves at the levels (/ + 1) and 
(m + 1). Notice that colored forests with (2rfc) white leaves, no black leaves, 
and coalescent degree r^, with different pairs of coalescent edges at level 
k, are necessarily of the form 

f(tfe) TT 7ptk,l.m 

k — 11 fc,/,m 
k<l<m<n 

for some families of integers = {tk,i,m)k<i<m<n such that ||tfc|| =Y.k,i,mtk,i,m = 
Tfc. Therefore, a colored forest in J^-q, with coalescence degree ||q||/2, a co- 
alescence sequence r such that ||r|| = ||q||/2, without a trivial tree with a 
white leaf, has necessarily the following form: 

f (t) _ f(to)r7-ro f(*")r/'»'n 

for some sequence of families of integers t = (tfc)o<fc<n5 such that ||tfc|| = r^, 
for any < /c < n. In the above displayed formula, Uk denotes the unique 
trivial tree with a single black leaf at level k, with < A; < n. We write t! 

fol' ]\k<l<m<n'^k,l,ra-- 

We are now ready to extend the Laurent expansions, and the Wick formula 
presented in Theorem 3.14, to particle models in path-spaces. 

Theorem 4.9. For any q= (go, ••• with ||q|| < N and any F e 
;B^^™(E^), we have the Laurent expansion 

(4.2) ^1^{F) = L'^+ E ^u^'^n^F) 

\ l<fe<||(q'-l)+|l / 

with the signed measures S'^'Q^ q defined by 

r<qM|r||=fcfgjF_(r) '^Ifl 

In addition, for any even integer ||q|| < N , and any symmetric function 
FGi3o^°'(E^), we have 

V/c < ||q||/2 d%,,^^{F) = and 

r<qM|r|| = ||q||/2(t>=r 

with the integer sequence (t) := (||tfc||)o<fc<n) '^i^'d the sum of the diagonal 
terms 6{t) = J2o<k<i<n'l'k,i,i- For odd integers ||q|| < A'', the partial measure- 
valued derivatives are the null measure on B^'^{'E^), up to any order 

k<m\/2\. 
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Before getting into the details of the proof, we mention that the Wick 
formula stated above has a natural interpretation in terms of the Gaussian 
fields (Vfc)o<fc<n introduced in Section 3.6. Following the discussion given in 
that section, we get that 

for any tensor product function F of the following form: F = -Fo.go ^ ' ' " ® 

with = ^ E.,6S,, (V'?^'^® • • •®</'fc^^"=^), G B,{Eu) and 7fc(v^D 
for any I <i < qk and any <k <n. 

Proof of Theorem 4.9. By definition of the particle model, and ar- 
guing as in the proof of (2.1), we find that 

IE([(7n^_i)^''"-®(7^)^''1Miei-l) 

qn 

for any ip G ""(^^r/ x E^"). This yields that (recalling Definition 4.7) 

w^(i^)iea)=7„^L^Q:iD^,. (F), 

from which we readily conclude that 

E(7;?''^(i^))=IE(7^^[Q^D^^, (F)]). 

A simple induction yields that 
E(7^'^(i^))=]E(^o®'°^'°D^^, Q?D;>^,...Q^D^^, (F)) 

iViq'l 4± (q')lal '"^ NWA^ (q')i?i ''^ 

from which we find the following formula: 



l<p<q' ^^^IpI fe:^_;|f|=p 

Using the Stirling formula (2.2), we readily check that 



^,.= E E ^(P'i)]^7^ E #(fXq 

l<l<q'l<p<q' ^"I^IpI fe;^_;|f|=p 
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From previous computations, we conclude that 

C=E]^ E K|fU'-r)7-^#(f)Al,. 

Finally, we notice that J-"q(0) reduces to the single class of all sequences of 
bijections in ^q. The end of the proof of the first assertion is now clear. To 
end the proof of the theorem, we notice that [with t as in the expansion of 

aii<iii/2Q„„(F)| 



<l<m<n 

Since we have s(|f^*^|,q' - r) = 1, and (q')|f(t)| = I\-i<k<n{<l'k)q'^-rk+i^ 
any r = {rk)o<k<n-, we conclude that 

s(|f(*)|,q'-r)^^-(t),_ q! 



(q')|f(t)| ' 2^(t)t! 

with t! = Y{Q<k<i<m<n^k,i,rn-- The end of the proof of the theorem is now 
straightforward. □ 

4.6. Propagations of chaos-type expansions. This section is essentially 
concerned with applications of the differential forest expansion machinery 
developed earlier to propagations of chaos properties of interacting particle 
models. In order to state and prove the main results of this section, we need 
to introduce some notation. We shall work throughout with a fixed time 
horizon n > 0, and a constant particle block size q < N. We let N^"*"^ be 
the set of integer sequences p = {pk)o<k<n £ N""*"^ such that ||p|| = q. We 
associate to a given p = {pk)o<k<n £ N""*"^ the pair of integer sequences p' 
and p + q defined by 

yO<k<n p'j^:= ^ pi and p + q := {po, ■ ■ ■ ,Pn-i,Pn + q)- 

k<l<n 

We also denote by {Gk)o<k<n and the collection of functions defined 
by 

Gk:=^^,x{r,k{Gk)-Gk) and := ® • • • ® Cf" 

(recall the G^'s have been defined in Section 1). Finally, we consider the 

integral operators Q,^^'^''^ from E^"'''' into E^_^i, defined for any function 
Fn+i G Bi){E'^_^j^) by the following formula 



44 

with 



P. DEL MORAL, F. PATRAS AND S. RUBENTHALER 



F ■=-^ 

7^(1)" 



(F„-,?^'?(F„)). 



The fohowing proposition and lemma are technical results that will be 
needed in the proof of Theorem 4.12. 

Proposition 4.10. For any q< N and any n>0, we have the polyno- 
mial decompositions 



:= E((l - 7^(G„)/7n = 



q/2<k<{n+l)(q-l) 

The derivatives of order q/2 < k < {n + l){q — 1) are given by the following 
formula: 

pgNj+ifc<||(p'-l) + |j P- 

Proof. We first check the following decomposition: 
(4.3) l-7^(G.)/7n(G„)= E ^p(Gp). 

0<p<n 

To prove this formula, we notice that 

7;^(G'„)-7n(G„)= n VpiGp)- n MGp) 

0<p<n 0<p<n 

+ [7^„i(G„_i) - 7n,-i(G'„_i)] X r?„(G„) = • • • 
= E lpiGp-Vp{Gp))x n VkiGk] 

0<p<n Lp+l<k<n 

with the convention Yl^s = 1- Finally, combining the multinomial decompo- 
sition 



(4.4) 



(l-7^(G'„,)/7n(G„))'^= E ^yT^'^lG^) 



with the Wick expansion stated in Theorem 4.9 we conclude that 



E 



q,n 



^ p! 

1^ peN"+i9/2<fe<||(p'-l)- 



E 



peN" 



^^d'^QnJG^n)- 



This completes the proof of the proposition. □ 
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Lemma 4.11. For any u G R — {1}, and m>l, we have the decomposi- 
tion 

^ 0<k<m l<k<q+l ^ / V\i U/ / 

with {q + k)k = {q + k)\/q\. 

Proof. This formula seems classical but we did not find it in the lit- 
erature. The proof is essentially based on the fact that for any n > and 
ti 7^ 1 , we have 

1 „,n+l pin f 

This implies that for any m> n, we have 

_ n! _ ^ fu-^^^ 
" (1 - n)"+i " ^du-\l-u^ 

^ ' n<k<m 

Applying the Leibniz binomial derivation formula for -§^{fg) to the product 

m.+l X _!_ ^i^g,^ 
l—u 

The lemma follows. □ 



Theorem 4.12. For any n>0, the sequence of probability measures 
(JPn+i,q)N>q is diffcrentiablc up to any order with d^¥n+i^q = iq^l^, and the 
partial derivatives given by the following formula: 

ofcp _ V- ((g-l) + |p|)! fc;^ -^(p+g) 

P^U()<i<2fc"i 

At any order N >q, we have the exact formula 

l<k<l{N-q)/2] 

with a remainder measure I^^+i^g such that sup^>]^ A^((^~'^^+^)/^||M^_^;^^^||tv < 
oo. 



Remark 4.13. Contrary to the expression in Theorem 1.5, where we 
have an exact expansion as a sum of polynomial terms, the expansion above. 



46 P. DEL MORAL, F. PATRAS AND S. RUBENTHALER 

though exact, is the sum of a polynomial term plus a remainder term which 
is not polynomial. This particular form comes from (4.5) appearing below 
in the proof of the theorem, where we have to develop a fraction, whose 
expansion is an infinite sum. 

Proof of Theorem 4.12. We let Fn+i be a bounded measurable sym- 
metric function on E^^;^, such that r7®^;^(F„+i) = 0. By definition of the 
particle model, we have that 

with the sequence of random variables := (1 — 7^ {Gn) /^n{Gn))- Using 
Lemma 4.11, we find that 

0<fc<m ^ InK'^n) / 

with the remainder term 

We also have from Proposition 4.10 

supjv>i ^/]VE(|7^(G„) - ^uiGnT^^f'"^^^ < 00. 

Using the regularity hypothesis (1.1) on the potential functions we conclude 
that 

SUp^.>i7V(-+l)/2p^^||TV<00. 

Finally, using the multinomial decomposition (4.4), we conclude that 

0<fc<m ^ ln{^n) J 

0<k<m pgf^^+i 

This yields, for any m + q < N , the functional expansion 



■ n+l,q fo— 1)!d! ^",P+<?^" "t -'^m 

||p||<r?i ^ ' 
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with a remainder measure -R^^, such that supjv>]^ A^*^™'"*"^^/^||i2^^||TV < oo. 
This imphes that for any k < (m + l)/2 we have 

||p||<m ^'^ 

This proves the first assertion of the theorem. Notice that feth-order deriva- 
tive measure d^Q^^p^q only involves colored forests with less than k coales- 
cent branches, from the original root, up to the final level. If ||p|| > {2k), 
then these colored forests contain at least one elementary tree with a white 

leaf. By definition of the operator Q^^^'^'' we find that 

This yields the second part of the theorem. □ 

5. Consequences. We end this paper with a series of some direct con- 
sequences of the above theorem. Some of the already known properties of 
particle systems can be found in Section 1. 

• The differential forest expansions presented in Theorem 4.12 allow us to 
deduce precise strong propagation of chaos estimates. For instance, for 
any N>{q + 7) [so that { [{N - q)/2\ - 1) > 2] we have 

1 

N' 



sup A^^ 



71+1,1? '/n.+l j^,-'-' '^n+1,. 



< oo 

TV 



N>q+7 

with a first-order partial derivative given by the formula 

d'F^+l,q{Fn+l)=d'Qn,qQT+l(Fn+l)+q ^ 5^Qn,q™ (i^n+l) 

0<m<n 

with the sequence of integers = (Im(fc) + q^n{k))o<k<n- The first term 
in the right-hand side in the above displayed formula has been treated in 
Corollary 3.13, and we have that 

d'Qn,qQnll(Fn+l) 

q{q - 1) ^ A fi,fe^®g 

- 5 ^n,q(4n+l\^ n+l) 

0<k<n 

0<k<n •'^k 



Each of the terms 5^Q„^q^Q^'"(F„+i) only involves the colored forests 

VO < A; < m lk,m ■= UkTk,m,nU^+l 



48 P. DEL MORAL, F. PATRAS AND S. RUBENTHALER 

associated with the trees Tk^m,n, and Ui, introduced on pages 816-817. 
After some elementary manipulations, we find that 

9^Qn,q™Qr(Fn+l) 



= q ASCQ^'"(i^n+i) 

0<fc<m 

0<k<m -^fe 

xQS+i(F„h.i)(x\...,x''). 

• A Wick formula derives from Theorem 4.9. More precisely, for any F £ 
B^"^ {E'^_^_^) , the partial derivatives d^Fn+i^q{F) are null up to any order 
k < [q'/2J , and for any even integer q, we have 

(5.1) 9''/¥„+i,,(F) = 7„+i(1)-''9''/2q„,,Q^|,(F) 

as soon as > 2{q + 2) [so that {[{N - q)/2\ - 1) > q/2]. To prove this 
claim, we notice that for any k < q/2, and any F G Bq^"\E^^i), we have 

VpG U Nf+\ /c<g/2<(||p||+g)/2 and therefore 

0<l<2k 

This yields that d^¥n+i.q{F) = 0, for any even integer g, and any k < q/2. 
In the case k = q/2, we have k = q/2 = (||p|| + q)/2 if, and only if, p 
coincides with the null sequence 0. 

• We let ll/illgsym = supj^g^sym |//(F)| be the Zolotarev seminorm on Ai{E'^^^) 
associated with the collection of functions 

B^^r--={FeBr{El+i)--\\F\\<l}- 
For any even integer q such that {q + 2) < N/2, we have 

,5.2) ^™p^A'-./a p. _ _i__l_a,/.Q„,,Q«., 



< oo. 



„sym 

"o,i 



• Combining the Wick formula stated above with the Borel-CanteUi lemma, 
we obtain for all q > 4 the almost sure convergence result 

(5.3) hm (7?^)®«(F)=0 p.s. 

N^oo 

for any bounded symmetric function F S Bq^^{E'/^). This result is an ex- 
tension of the law of large numbers for [/-statistics obtained by Hoeffding 
[8] for independent and identically distributed random variables to in- 
teracting particle models. One can also look at [9] for a more modern 
reference. 
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We mention that the same hnes of arguments used in the proof of Theo- 
rem 4.12 show that the sequence of probabihty measures 



is differentiable up to order [{N — q)/2\ , with d^Fn+i,q = Vn+i^ 

partial derivatives given for any 1 < k < [{N — q)/2\ by the fohowing 

formula: 

In the same way, for any F £ BQ^^{E!i[^^) and N > 2{q + 2), the partial 
derivatives (9^Q„+i^q(F) are null up to any order k < q/2, and for any 
even integer q, we have that 

d''/^F„+i,,{F) = 7n+i(l)-''a'?/2Q„+i,,(F). 

• Finally, the Wick formula stated above allows to deduce sharp Lg-mean 
error bound. To see this claim, we simply observe that 

F = (/-r?„+i(/)r 
with / e =^ P^+i,,(F) = E([r?^^+i(/) - i]n+iim 

and for any even integer q such that (g + 2) < N/2, we have the inequality 

(5.4) sup N^+^Z^ f^,^^ ^— L— S-'/^Q , 



N>q+7 



< oo. 



„sym 
"0,1 
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